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Abstract 

An eikonal model has been used to assess the relationship between calculated 

strengths for first forbidden /3 decay and calculated cross sections for (p, n) 

charge exchange reactions. It is found that these are proportional for strong 

transitions, suggesting that hadronic charge exchange reactions may be useful 

in determining the spin-dipole matrix elements for astrophysically interesting 

leptonic processes. 
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I. INTRODUCTION 



Charge exchange reactions {Az, Az±i) induced by hadronic projectiles are a powerful 
tool for probing spin-isospin degrees of freedom in nuclei The spin-isospin parts of 

the operators that mediate charge exchange reactions such as {p, n) are the same as those 
involved in the corresponding processes induced by electromagnetic and weak interactions. 
As a result, the matrix elements that describe hadronic charge exchange reactions are closely 
related to those that describe the rates of /3 decay or the cross sections of reactions induced by 
neutrinos. It would be fortunate if this relationship were quantitatively accurate, since it is 
often difficult to study the leptonic processes directly. For example, the range of excitation 
energy kinematically accessible in a /? decay transition does not encompass the majority 
of the allowed (Gamow- Teller) strength and the experimental study of neutrino induced 
reactions is difficult. 

A promising direction of future activity is to determine leptonic strengths for otherwise 
inaccessible nuclides by studying charge exchange reactions using radioactive (secondary) 
beams in inverse kinematics 0. This would provide nuclear properties important for prob- 
lems of nuclear physics, particle physics, astrophysics and cosmology One could 
clarify the relationship between the spatial properties of nuclear halo systems and the na- 
ture of soft multipole modes |^ . One could also determine the strength of neutrino-nucleus 
interactions needed to describe the chemical evolution of the Universe, especially the abun- 
dances of the light elements |^ and the products of r-process nucleosynthesis 0], and to 



calibrate terrestrial detectors of supernova neutrinos |[Tl 



However, it is not obvious a priori that the correlation of charge exchange and leptonic 
matrix elements is sufficiently close for this purpose. In contrast to leptonic processes, 
hadronic reactions involve operators that have an additional radial dependence and are sub- 
ject to distortion by complex nuclear potentials; the medium renormalization of effective 
operators and the contributions of multi-step processes introduce additional uncertainties. 
Therefore it was an important advance to establish that there is an approximate proportion- 
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ality between the cross section of charge exchange reactions at very forward angles leading 
to the Gamow- Teller (GT) excitations (with transferred T = 1, L = 0, J = S = 1) and 
the transition strength B{GT) determined by intrinsic nuclear matrix elements PJT^-[I6 



The proportionality has been confirmed for strong transitions in a variety of nucleon and 
nucleus induced charge exchange reactions; a more detailed analysis is needed for weak GT 



processes [|T^. The model-independent character of the relation between the L = cross 
section for reactions induced by ^^C projectiles and the GT strength was clarified by a the- 
oretical analysis based on a sensitivity function which identified the important part of 
the target transition density in momentum space. 

In contrast to GT transitions, the first forbidden matrix elements of weak processes 
explicitly include orbital degrees of freedom. The corresponding nuclear response in the 
L = 1 channel is associated with the states forming the spin-dipole and giant-dipole charge- 
exchange resonances (SDR and GDR). This excitation was discovered and studied on 
different targets |^ mostly using {p, n) reactions. Recently the energy splitting of the 
L = 1 charge exchange resonances (GDR and SDR) was determined [^. There is very 
little information about the quantitative relationship of charge exchange cross sections and 
leptonic strength for these excitations. Here we take a first step in providing this information 
by studying the relationship between first forbidden strength and {p, n) cross sections, both 
calculated from the same wave functions. We show that, at the same level of accuracy as for 
the GT case, for strong transitions one can expect an approximate proportionality between 
the observed cross sections of charge exchange reactions populating spin-dipole states and 
the corresponding nuclear transition probabilities. We follow the general approach that was 
successfully applied to GT excitations by Osterfeld et al. [|18[, extended to describe L = 1 
transitions and the effects of the real part of the optical potential. 
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II. THEORETICAL FRAMEWORK 



In order to understand the relationship between the charge exchange cross section and 
the nuclear response strength, we need to examine the effects specific to the excitation of 
the SDR. For this purpose we consider the influence of distortion in a simple eikonal approx- 
imation (EA). In many cases, even at rather low energy, the EA gives a good qualitative 
description of the reaction cross section in the SDR region. We obtain the relevant wave 
functions and transition densities from the shell model. 

The strength Bj of the SDR in the long wavelength limit and the transition form-factors 
Fj{q) are calculated in terms of the elements of the single-particle density matrix of the 
given transition i ^ f from the ground state, p/j(z/, u') = {f\al,a^\i), 



Bj{i - /) 



^pfi{v, p'){p'\\rOj\\v) 



(1) 



Fj{q) = Y.Pf^{^^^')Wmqr)O.J\\u), (2) 

vv' 

where {i''\\ji{qr)Oj\\v) is the product of the radial matrix element of the spherical Bessel 
function and the reduced matrix element of the charge exchange spin-dipole operator Ojm = 
{cr (g) Yi(f)}jjv/r^. In the limit of low momentum transfer, qR <^ 1, the squared transition 
form-factor is related directly to the strength of the SDR, 

FKq)^JBJ{^^f). (3) 

In the case of the GT resonance a direct proportionality between the experimentally 
measured cross sections of charge exchange reactions at very forward angles and nuclear 
matrix elements was confirmed by a number of studies p|,|l2|-^] at a level of accuracy of 



10-15%. From the viewpoint of the underlying physics, this important result is based on a 
single-step mechanism for the process, a simple bare operator which does not include orbital 
degrees of freedom, and the dominance of the central spin-isospin interaction Var over a 
broad range of energies. It is a priori unclear whether these features pertain to the SDR 
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case (AL = 1). At the maximum of the differential cross section for the SDR, qR ~ 1. The 
g-dependence of the hadronic operator and the effects of tensor forces may show up at the 
larger q. As a result, there may not be a simple relationship between the cross section and 
the nuclear response strength. 

In the SDR case, the direct part of the (p, n) reaction amplitude (the exchange part will 
be discussed later) is 

TfT = J dhx^->{kf,r) J ^F,(g')K7M(q')exp(-.q'-r)xS+)(k„r), (4) 

where the effective operators for the channels with angular momenta J'^ = 0^ , 1^ and 2~ 
contain contributions from central and tensor forces, 

VjMic{) = Vj\,ic{) + VjM), (5) 



VjM) = Am^ ^^tliq')[{a ® r;(q')}jM - 3(^t ■ q'){q' ® Y^^jm]. (7) 

In Eqs. (|]) and (^, f^^-iq') and t* (g') are, respectively, the central and tensor components 
of the nucleon-nucleon t-matrix, see Franey and Love that are responsible for spin- 
isospin transfer. The excitation of different J-components of the SDR proceeds via different 
combinations of the amplitudes of the nucleon-nucleon effective interaction. For the 0~ part, 
the tensor interaction can be combined with the central one using the relations 

{q'®n*(q')}oo = --^, {a®Y*moo = -^-^jfi. (8) 
v47r v47r 



These identities produce the combination = t^^-iq) — 2t* (g) which is just the spin- 
longitudinal component of the nucleon-nucleon t-matrix. For the 1~ part, the operator {q'® 
Y{{q')}iM = 0, and the amplitude in Eq. (|) becomes proportional to = t^^(g) + t* (g), 
which is the spin-transverse component of the nucleon-nucleon t-matrix. For the 2^ part, 
both components contribute, and the amplitude in Eq. has a more complicated form. 



The functions x^~^ X^^^ the amplitude (^) are optical-model wave functions de- 
scribing the motion of the initial proton and final neutron in the optical potential of the 
target nucleus. To disentangle the nuclear transition form-factor from the observed cross 
section, one needs to unravel the intrinsic nuclear dynamics masked by the distorted waves 



III. DISTORTION FACTOR 

The effective operators and (|^ in the reaction amplitude are evaluated at the 
value q' of the local momentum transfer that corresponds to the charge exchange event. 
However, because of the distortion by the optical potential, q' does not coincide with the 
asymptotic momentum transfer q = kj — kj. In the absence of distortion (the plane wave 
approximation) we would have 

[X/"^*(k/, r)xS^^(ki, y)]pw = exp (zq ■ r) (9) 

so that the integration over r could be performed explicitly resulting in 5(q — q'). For 
distorted waves this is no longer true. 

In the EA, the product of two optical-model wave functions in Eq. @ can be estimated 

by 

X^"^*(k;,r)xS+^(k„r) = exp(zq-r)D(r^), (10) 



where the distortion factor D{r±) is defined by 



D{r±) = exp 



T / —Uopt{z,Yi_) 
Tl J-oo V 



(11) 



In the spirit of the eikonal approximation, the longitudinal momentum is still preserved, 
while the distortion is effective in the plane perpendicular to the trajectory. In Eq. (pTD, 
the optical potential is different in the initial and final channels. In first order we account 
for this difference by assuming a fast single-step process which leads to 
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TT 1 //"P TT^ ' 

opt J- / ^ opt opt 



(12) 



V 2 y t>j t>/ 

For a square well potential of depth Uq with a sharp boundary at r = i? the distortion factor 
can be calculated analytically: 

2f/f 



D(r) = exp ^ - r2 , (13) 

\ hv ) 

and D(r) = 1 for r > R. This approximation might be insufficient in the region of minimum 
of the cross section where the details of the potential shape are essential. However, near the 
maximum, which is our region of interest, the cross section is insensitive to the diffuseness 
of the optical-model potential. We have checked this point by varying the diffuseness in a 
DWIA calculation. 

The exchange part of the reaction amplitude was estimated in the "fixed Q" approxima- 
tion [^,^. This approximation has been found |2|j2|] to improve with increasing energy 



and multipolarity of the excitation for even-state forces and mainly S* = natural parity 
transitions. We need to deal with S = 1 transitions and include tensor forces. Recently, the 



fixed Q approximation for exchange was discussed in detail in a review article ||2^ ; it is more 
accurate for central interactions than for the tensor interaction so that the latter is of most 
concern. Fortunately, the approximation for tensor forces is best when Q is not far from /cj, 
which is the case for the relatively high energies considered here (fcj is greater than the Fermi 
momentum in the nucleus) and for small angles. It is also relatively less important for the 
transitions we consider where there are large central and direct tensor contributions ||28|| . 
Baker et al. |^ note that for the intermediate energy {p,p') isovector excitation of 1^ and 
2~ states in ^°Ca, the cross sections obtained with the fixed Q approximation were usually 
within 30% of the results of calculations done with exchange treated exactly. As shown in 
Fig. 1, our eikonal calculations with the fixed Q approximation are within 20% of DWBA 
calculations with exact exchange, consistent with the {p,p') result. This level of agreement 
seems adequate for our purposes: the overall effect of the approximation is simply to change 
slightly the scale of the cross sections for a given J from that shown in Fig. 2. If the effect 



were different for different transitions, it would show up in the overall scatter of the points 
about the average line which is small for strong transitions. 

In the fixed Q approximation the exchange momentum Q coincides in the laboratory 
frame with the initial momentum kj. For the exchange amplitude we then obtain 



(14) 



where the central t^^, and tensor, t^, interactions are defined by Franey and Love [^. In 
Eq. (|l^ the effective exchange operator 



(15) 



includes the distortion factors D{r) specific for each nucleon inside the nuclear matrix ele- 
ment. 



IV. SENSITIVITY FUNCTION 

In Eq. (^) the integration over r is not well defined at large distances. It is convenient to 
single out a no-distortion contribution proportional to 5(q — q') by using the decomposition 
D{r) ^ 1 + [D{r) — 1]. The first term describes the plane wave contribution and the second 
the effects of distortion. Since |-D(r)| < 1, the distortion term reduces the plane wave 
contribution. A convenient form of Eq. (^) can be obtained for transitions to 0^ and 1~ 
states by writing it as 



~ V2J+1 ' ^ ^ 

where 

Ti^^ = TiiPW) + / dq' Siiq, q')Fj{q') (17) 
Jo 

is the amplitude describing the excitation of the SDR with the longitudinal, m = 0, or 
transverse, m = ±1, relative proton-neutron spatial oscillations. In Eq. ([T^), T^\PW) is 
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the plane wave contribution, and we have introduced the sensitivity function |T8 

Si{q. q') = f dh exp (zq ■ r){D{r^) - l)Y,Ur)jiiq'r) 

I t'(g') for J = 0, 
X . (18) 

[ t'^{q') for J = 1 

The sensitivity function in Eq. (0) characterizes the range of q' which contribute impor- 
tantly to the charge exchange cross section for a given asymptotic momentum transfer q. 



V. EXAMPLE: ^"^CiP, N)^'^'N REACTION 

As an example of application of the method we performed numerical calculations for the 
^^C(p, n)^^N reaction to compare with experimental data for the excitation of spin-dipole 
states at a proton energy of 135 MeV. 

For ^^C, with the optical- model potential of Ref. |^0|, D{r) varies smoothly inside the 



nucleus. Near the surface it changes rapidly from its value at the center -D(O) ~ 0.5, to the 
value of 1. It is then a good approximation to write the exchange matrix element in Eq. 



(0) as 

(/|0'^%)|z) ^D(ro)(/|Eexp(zq-r,Kr-|z). (19) 

3 

The result is not very sensitive to a particular choice of the reference point ro; we used 
tq = 0. This approximation underestimates the exchange part. Near the maximum of 
the cross section its contribution is not significant. It becomes important at large angles 
where the difference of distortion along different trajectories is noticeable; however, the cross 
section at large angles is small. We mention parenthetically that, opposite to the results 



of p5|j26[1 for even-state forces, the exchange contribution reduces the cross section. This 
is related to our inclusion of odd-state forces whose net effect is small after exchange is 
accounted for. 

In our calculations, the wave functions and transition densities for the spin-dipole states 
were obtained using a harmonic oscillator basis including the orbitals of p, sd and pf shells 



that form the 3 hoo model space necessary for the description of the L = 1 excitations. The 
calculations were performed with the WBN residual interaction and a harmonic oscillator 
parameter of 1.64 (A/A - 1)^/2 fm [0. The cross sections for the ^^C(p, n)^^N reaction 
leading to the 1~ state at = 1.8 MeV and the 2~ state at = 4.3 MeV were calculated 
as the sum of the direct and exchange amplitudes, Eqs. (|^ and (|I^). The results are shown 
in Fig. 1 together with data from Ref. ||29|. For comparison, a calculation with the DW81 
code is also presented. The calculations give similar cross section shapes near the maximum 
and significantly overestimate the magnitude of the cross section. The results are very 
similar for other excited states. They are also similar to the distorted wave results obtained 
for the same transition in Ref. ||29| using a Ihu model space and the MK interaction. 

The systematics of the cross sections at their maximum divided by the calculated (3 decay 
strengths are shown in Fig. 2 for 0^, 1^ and 2^ states. As seen from Fig. 2, there is an 
approximate proportionality between the cross section at the maximum and the spin-dipole 
strength, accurate to within 10-15%, for states with strength Bj > 0.1 fm^. This is the 
same level of proportionality as for GT (L = 0) excitations at very forward angles. One 
may ask whether the validity of this conclusion is affected by the poor agreement in the 
magnitude of the cross section for the 1.8 MeV 1~ state. We would argue that this is not 
the case: since the wave functions are sufficiently complex, they provide a reasonable sample 
of possible behavior with respect to the operators involved. Furthermore, one might expect 
proportionality to fail for such weak transitions. 



VI. DISCUSSION 

It is not clear a priori that the high degree of proportionality shown in Fig. 2 should 
occur. The cross section involves an integral of the transition form factor over a range of g' 
while the value of Bj is determined by evaluating the form factor at q' ~ 0. To examine what 
leads to the observed proportionality, we return to Eq. ((T^). Two factors determine the 
(p, n) cross section: the transition form-factor Fj{q) and the sensitivity function S^{q,q'). 
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In Fig. 3 we show the transition form-factors for different 1~ states normahzed to the 
same maximum value in order to compare their shapes. The shapes are very similar near 
the maximum but differ at higher momentum transfer q. If the region of high q' does not 
contribute significantly in the integration over q' in Eq. (p!7D, the integrals for different 
form-factors will be proportional. 

Samples of the imaginary parts of the sensitivity functions are shown in Figs. 4 and 5; 
the real parts have very similar shapes and are typically a factor of two smaller in magnitude. 
A general remark should be made about the g'-dependence at small q' . Since D{v±) does not 
depend on the longitudinal coordinate the 2;-component of the local momentum transfer 
q' must coincide with the 2;-component of the asymptotic momentum transfer q. When the 
absolute value of q' is smaller than g^, this condition cannot be fulfilled at any angles of q', 
and the sensitivity function must be equal to zero. 

As noted above, the sensitivity function for J = is proportional to the spin-longitudinal 
component of the effective interaction, and that for J = 1 to the spin-transverse one. We, 
therefore, expect a different g'-dependence reflecting the different behavior of t\q') and 
t*^(g'). For 0~ states the projection m = dominates, corresponding to the spin- longitudinal 
behavior of the reaction amplitude for J = 0; the sensitivity function for m = 1 is smaller 
by an order of magnitude. At the small scattering angle corresponding to the peak of the 
cross section, 9 = 4.3°, the momentum transfer q is almost parallel to the initial proton 
momentum kj, thus enhancing the m = component. For J = 1 the picture is different, as 
is seen in Fig. 5. Projections m = and m = ±1 give comparable contributions. 

Given the nature of the sensitivity functions it is clear why the cross sections and Bj 
are closely proportional. For both J = and J = 1 the main contribution comes from the 
peak region where the transition form factors have the same shape, leading to the observed 
proportionality. 

In summary, our results imply that there will be an approximate proportionality of the 
observed cross section at the maximum of the charge exchange reaction exciting spin-dipole 
modes and the leptonic strength. This supports the possibihty of using such reactions for 
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extracting leptonic strengths of astrophysical interest. Having established here the basic ap- 
paratus to examine this issue, it will next be important to consider transition densities for 
heavier nuclides, so as to determine whether their shapes are similar enough that cross sec- 
tions and Bj strengths will be proportional. It will also be important to examine the nature 
of the sensitivity functions for heavier nuclei, to ascertain whether they remain concentrated 
in a relatively small range of q' where the transition densities are similar. 
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FIGURES 

FIG. 1. Cross sections for the reaction ^'^C{p,n)^'^N leading to the 1~ state at Ex = 1.8 MeV 
and the 2~ state at Ex = 4.3 MeV. The cross sections shown as solid lines are the results of 
the eikonal approximation calculations described here; a DWIA calculation done with the same 
parameters is shown by solid dots. The data shown as solid squares are from Ref. as are the 
DWIA calculations shown as open circles. All the theoretical calculations have been multiplied by 
the factor shown in the Figure. 

FIG. 2. Ratios of the cross sections for excitation of the spin-dipole states, taken at their 
maximums, to the corresponding spin-dipole strengths Bj for states with different Bj. The upper 
panel is for 0~ states; the middle panel for 1" states; and the lower panel for 2^ states. 

FIG. 3. Transition form-factors for the 1^ states normalized to their maximum values. The 
state with the anomalous shape corresponds to the high point near Bi = 0.45 fm^ in Fig. 2, middle 
panel. 

FIG. 4. The imaginary part of the sensitivity functions Sm{Q,Q') for J = 0. The small size of 
Si reflects the spin-longitudinal origin of the reaction amplitude. 

FIG. 5. The imaginary part of the sensitivity functions Sm{q,q') for J = 1. Si and Sq are 
comparable for this spin-transverse dominated reaction amplitude. 
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